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摘 要: 基于 CFD(Computational Fluid Dynamics ) 数 值 模拟 ,以 青海 省 中 灶 火 省 道 303 公 路 为 研究 区 域 , 对 流动 沙 
丘 区 公路 路 基 风 沙 流 场 进行 了 模拟 分 析 ,揭示 了 流动 沙丘 区 路 面 沙 害 的 形成 机 制 ,以 期 为 流动 沙丘 区 路 面 沙 害 的 
防治 提供 理论 依据 。 数 值 模拟 结果 表明 : 当 气 流 经 过 流动 沙丘 和 路 基 时 ,气流 速度 产生 分 区 ; 当 存 在 流动 沙丘 时 ， 
路 基 周 边 的 减速 区 范围 明显 增 大 ,路基 坡 面 和 背风 坡 全 部 处 于 气流 低速 区 的 范围 ,沙丘 高 度 以 上 ,在 沙丘 顶部 气流 
速度 急剧 增 大 , 受 障碍 物 自身 形状 的 影响 ,风速 变化 由 “0Q” 型 分 布 变 为 “M” 型 分 布 。 当 沙丘 表面 的 原 有 防护 体系 下 C 
于 维护 更新, 沙 粒 在 沙丘 附近 大 量 堆 积 ,沙丘 会 慢 慢 治 主 风 向 移动 ,成 为 沙 源 , 捧 埋 路 基 而 形成 沙 害 。 数 值 模拟 结 
果 与 现场 实际 相 吻 合 , 证 实 了 数值 模拟 的 准确 性 。 对 于 流动 沙 夺 区 公路 的 沙 害 防治 ,应 详细 考察 当地 的 地 貌 形态 
及 风沙 流 运动 规律 ,定期 检查 维护 原 有 防护 体系 ,及 时 清理 积 沙 , 既 要 阻止 风沙 流 对 沙丘 进行 沙 源 供给 ,又 要 防止 


沙丘 本 身 成 为 沙 源 。 
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风沙 活动 是 塑造 地 貌 形 态 的 基本 动力 和 物质 
基础 ,也 是 风沙 灾害 产生 的 主要 条 件 "。 只 要 风能 够 
使 地 表 沙 粒 运 动 和 堆积 ,就 会 对 承载 体 产生 沙 害 ” 。 
风沙 活动 不 仅 在 地 质 历 史 时 期 对 全 球 地 貌 形 态 的 
演变 起 了 重要 作用 ,而且 在 现代 文明 时 代 对 人 类 生 
存 环境 和 生产 活动 也 产生 着 巨大 影响 。 沙 漠 独 特 
的 气候 条 件 造 就 了 各 式 各 样 的 沙丘 , 且 沙 丘 表 面 的 
沙 物 质 极 易 被 风 吹 扬 , 成 为 风沙 灾害 的 主要 沙 
源 地 。 

沙漠 地 区 气候 干燥 ,降水 匮乏 ,植被 稀少 , 沙 源 
丰 亩 ,在 强劲 的 大 风 作 用 下 ,主要 形成 了 路 基 风 蚀 
和 沙 埋 形式 的 沙 害 ”” ,其 中 , 沙 埋 是 路 基 沙 害 的 主 
要 表现 形式 ,流沙 上 路 ,不 但 会 造成 路 面 磨 蚀 , 增 
大 养护 成 本 ;在 大 风 等 恶劣 天 气 ,还 会 产生 扬 沙 ,能 
见 度 大 大 下 降 , 极 易 引 发 交通 事故 ,甚至 会 阻 断 交 
通 。 此 外 ,流动 沙丘 不 断 侵占 公路 ,造成 道路 交通 
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中 断 ,大 规模 沙丘 压 埋 甚至 会 迫使 公路 改线 ,大 大 
增加 建设 成 本 。 

随 着 我 国 “ 西 部 大 开发 "和 “一 市 一 路 " 倡 仪 的 
陆续 实施 , 越 来 越 多 的 沙 区 公路 投入 运营 ,为 了 保证 
沙 区 公路 的 安全 运营 ,大 批 学 者 通过 风 洞 试验 "和 
野外 观测 ”的 方法 ,对 沙 区 公路 的 路 基 风 蚀 和 积 
沙特 征 进行 了 人 研究 。 鱼 燕 萍 等 在 风 洞 试验 中 选 
取 了 两 种 典型 路 基 , 通 过 分 析 不 同 坡 角 和 高 度 对 路 
基 流 场 的 影响 ,研究 了 路 基 断 面 的 风沙 输 移 -堆积 
过 程 ;Li 等 "在 风 洞 中 进行 了 缩 尺 试验 ,以 不 同 直 
径 的 沙 粒 和 不 同 坡度 的 路 基 为 模型 ,模拟 了 沙 粒 对 
沙漠 公路 造成 破坏 的 过 程 ; 韩 致 文 等 “ ”对 穿越 塔 
里 木 沙漠 的 公路 进行 研究 ,阐明 了 风沙 灾害 产生 的 
动力 条 件 、 物 质 基础 以 及 风沙 流 在 路 面 的 沉积 过 
程 。 但 是 风 洞 试验 和 野外 观测 耗费 大 量 人 力 、 物 
力 , 且 受 仪 名 设备、 客观 环境 的 影响 较 大 ,无 法 直接 
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达到 理想 条 件 。 随 着 计算 机 技术 的 飞速 发 展 ,CFD 
数值 模拟 "作为 一 种 新 兴 的 研究 手段 ,以 其 方便 
快捷 、 成 本 低 及 能 够 处 理 大 量 数据 的 优势 ,被 国内 
外 学 者 广泛 应 用 于 风沙 问题 的 研究 。 李 驰 等 “5 通 
过 数值 模拟 的 方法 ,人 研究 了 不 同 路 基 高 度 ,不同 边 坡 
坡 率 对 路 基 周 围 风速 流 场 分 布 的 影响 ; 石 龙 等 "" 基 
于 数值 模拟 ,揭示 了 路 基 周 于 风沙 两 相 流 的 运动 特 
性 和 路 基 的 积 沙 形态 。 谢 虎 雄 等 "应 用 数值 模拟 
手段 ,针对 不 同 沙 粒 粒 径 、 摩 阻 风 速 和 沙 粒 跃 移 等 
问题 进行 研究 ,揭示 了 沙漠 公路 两 侧 风 沙 流 场 的 空 
间 分 布 特性 。 武 生 智 等 采用 数值 模拟 的 方法 ,人 研 
究 了 沙漠 公路 路 基 的 横断 面 形 状 对 近 壁 流 场 的 影 
响 ; He 等 "通过 数值 模拟 的 手段 ,研究 了 公路 桥梁 
以 及 桥墩 附近 沙 粒 沉积 速率 的 变化 特征 。 

然而 ,以 上 研究 成 果 均 是 关于 风沙 流 对 公路 路 
基 产 生 的 风沙 危害 ,在 流动 沙丘 区 域 , 路 基 不 仅 受 
风沙 流 的 风蚀 和 积 沙 危 害 , 还 受到 流动 沙丘 迁移 的 
压 埋 危 害 , 但 是 有 关 流动 沙丘 区 公路 路 基 的 风沙 活 
动 却 鲜 有 报道 ,流动 沙丘 区 公路 沙 害 防治 措施 的 比 
选 尚 无 明确 的 理论 依据 下 。 因 此 ,本 文 以 青海 省 
乌 图 美 仁 乡 中 灶 火 东南 3~10 km 处 的 省 道 303 公 路 
为 研究 区 域 ,通过 数值 模拟 的 方法 ,分 析 流 动 沙丘 
区 公路 路 基 的 风沙 流 场 变 化 规律 ,以 期 揭示 流动 沙 
丘 区 公路 路 面 沙 害 的 形成 机 制 ,并 为 该 地 区 公路 风 
沙 灾 害 防 治 提供 理论 依据 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 

研究 区 段位 于 青海 省 乌 图 美 仁 乡 中 灶 火 东南 
3~10 km 处 的 省 道 303 公 路 ,研究 区 段 主 要 以 风 积 地 
貌 为 主 ,还 包括 零散 分 布 的 各 种 沙丘 ;气候 类 型 为 
高 原 大 陆 型 干旱 性 气候 ,降水 稀少 ,蒸发 量 大 , 量 地 
下 水 位 深 ; 植 被 稀少 ,零星 分 布 有 沙 拐 束 、 树 柳 等 天 
然 植被 ;风力 强劲 ,年 平均 风速 为 5~6 ms ,最 大 风 
速 超 过 30 ms'3, 主 导 风 向 和 主要 起 沙 风向 均 为 
偏 西风 ; 沙 的 来 源 大 多 为 新 月 形 流 动 性 沙丘 , 沙 源 
充足 。 研 究 区 段 公 路 两 侧 分 布 大 量 沙 丘 ,前 期 为 了 
防止 沙丘 前 移 压 埋 路 基 , 工 作 人 员 在 沙丘 表面 采取 
了 沙 障 固定 的 措施 ,但 随 着 时 间 的 推移 ,沙丘 表面 
HI Ba SPAS ZA itt PAE DP .更 新 ,迎风 坡 沙 障 长 期 遭受 
严重 的 风蚀 作用 ,背风 坡 沙 粒 大 量 堆 积 , 使 得 沙丘 
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又 有 了 流动 的 趋势 ,从 而 对 路 基 造 成 严重 的 沙 埋 危 
害 , 如 图 1 所 示 。 


图 1 人 研究 区 公路 路 基 沙 埋 示意 图 
Fig. 1 Sand burial diagram of highway subgrade in 


the study area 


1.2 计算 域 及 网 格 划 分 

考虑 到 沙丘 形状 对 模拟 结果 的 影响 以 及 建 模 
的 难 易 程度 ,将 流动 沙丘 简化 为 图 2 所 示 。 经 现场 
实际 调研 ,流动 沙丘 及 公路 路 基 的 基本 参数 为 : 沙 
丘 高 Sm, 长 28m, 路 基 高 2m, 宽 12 m, 且 左右 对 
称 。 流 动 沙丘 设置 在 距离 人 口 40 m 处 ,模型 利用 
Gambit 前 处 理 软件 进行 网 格 单元 划分 ,在 底部 壁面 
设置 10 层 边界 层 , 其 中 第 一 层 边 界 层 的 高 度 设置 为 
y+, 增 长 速率 为 1.2, 最 终 整 个 计算 域 网 格 数目 为 
200x10' 个 。 为 了 达到 一 个 相对 可 靠 的 结果 ,网 格 划 
分 时 ,底层 网 格 一 般 布 置 到 对 数 分 布 规律 成 立 的 范 
EIN ,这 就 要 求 11.5~30<y+<200~400。 本 文 模拟 
中 ,y+ 的 最 大 值 为 110。 对 于 边界 条 件 的 处 理 , 采 
用 速度 入 口 (Velocity-inlet) 、 出 流 条 件 (Out-flow) 的 
方法 ,流动 沙丘 和 路 基 采 用 wall 边界 条 件 行进 描 
述 ” 。 入 口 处 风速 廓 线 取 对 数 分 布 方程 ,轴线 速度 
分 别 为 12 ms 18m.s 24m:s 和 30m:s 。 
1.3 研究 方法 

采用 欧 拉 双 流体 模型 对 风沙 两 相 流 运动 进行 
瞬 态 模拟 计算 ,将 风沙 两 相 流 运动 中 的 每 一 相 都 视 
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图 2 流动 沙丘 与 路 基 人 简化 模型 
Fig.2 Simplified model of mobile dune and roadbed 
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为 填充 在 整个 计算 区 域 的 连续 介质 ,在 同一 时 间 和 
同一 空间 内 进行 复杂 的 耦合 运动 。 流 场 模 拟 过 程 
中 忽略 温度 的 影响 ,将 流动 视 为 绝热 流动 ,同时 ,由 
于 进口 处 马赫 数 较 低 (来 流速 度 最 大 设 定 30 m+s”), 
认为 空气 压缩 性 对 气流 流动 特性 无 影响 ,为 不 可 压 
缩 气 流 。 所 在 整个 计算 域 为 求解 域 ,在 整个 求解 域 
内 满足 质量 动量 以 及 能 量 守恒 。 考 虑 到 K-epsilon 
模型 是 在 高 雷诺 数 条 件 下 发 展 的 计算 模型 ,而 在 近 
壁面 区 域 ,雷诺 数 较 低 , 故 在 近 壁 面 区 域 采 用 标准 
壁面 图 数 (Standard Wall Functions) 来 保证 计算 精 
度 ,采用 剪 切 压力 传输 模型 (SST) 开 展 流 场 计算 , 采 
用 二 阶 隐 式 方法 进行 时 间 积 分 ,采用 Phase Coupled 
SIMPLE 算法 进行 求解 ,清流 动能 .动量 等 离散 化 问 
题 均 采用 二 阶 迎 风格 式 。 计 算 过 程 中 残 差 控制 在 
10° ,达到 收敛 ,可 保证 模拟 结果 的 准确 性 。 
1.4 计算 参数 

经 现场 实际 调查 , 积 沙 主要 以 细 沙 和 极 细 沙 为 
主 , 故 在 Fluent 中 设 定 沙 粒 粒 径 d=10” m, 沙 相 密 度 
p:=2650 kg .m ,空气 密度 p=1.225 kg.m ,黏度 人 = 
1.79X10° Pa.s ,压力 设 定 成 常 压 ,又 因为 多 相 流 理 
论 为 稀 相 , 故 沙 粒 体积 分 数 设 置 为 2%。 动 力 粘度 、 
磁 撞 粘度 均 选 用 Gidaspow 模型 ,摩擦 精度 选用 


(a) V=12 m-s” 
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Schaeffer, 沙 床 最 大 堆积 率 为 0.63 ,并 考虑 相间 作用 
力 ,升力 系数 和 阻力 系数 均 为 saffman-mei, fil fe PK 
复 系数 为 0.9 ,动量 交换 系数 采用 Syamlal-0 "Brien 模 
型 。 


2 结果 与 分 析 


2.1 数值 模拟 风沙 流 场 分 布 特征 

由 不 同 风速 下 流动 沙丘 与 路 基 周 围 流 场 速 度 
可 知 (图 3), 不 同 风 速 作用 下 ,流动 沙丘 与 路 基 周 围 
风沙 流 场 分 布 特征 基本 一 致 。 风 速 V=24 ms THY, 
风沙 流 运动 到 沙 顾 附近 ,受到 沙丘 的 阻 涡 作 用, 风 
速 逐渐 降低 ,在 沙丘 迎风 坡 形成 气流 减速 区 , 且 速 
度 降幅 较 大 。 随 后 气流 沿 坡 面 候 升 ,在 沙丘 顶部 不 
汤 汇 集 , 形 成 气流 加 速 区 。 此 外 ,沙丘 迎风 坡 坡 面 
风速 等 值 线 与 沙丘 迎风 坡 坡 度 基本 平行 。 气 流 继 
续 向 前 运动 ,在 沙丘 背风 坡 .路 基 顶 面 及 背风 坡 附 
近 形 成 气流 低速 区 , 当 气 流 完全 通过 障 但 物 , 则 慢 
慢 恢 复 至 来 流风 速 。 

不 同 风 速 下 路 基 周 围 流 场 速 度 分 布 如 图 4 所 
示 。 对 比 图 3 和 图 4 可 以 看 出 ,流动 沙丘 的 存在 对 
路 基 周 边 的 流 场 分 布 产 生 了 显著 的 影响 , 当 流 动 沙 
丘 不 存在 时 ,气流 加 速 区 处 于 路 肩 处 ,路基 坡 面 处 


(b) V=24 m-s 
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图 3 不 同 风 速 下 流动 沙丘 与 路 基 周 围 模 拟 风沙 流 场 分 布 


Fig. 3 Distribution diagram of simulated aeolian sand flow field around mobile dune and roadbed at different wind speeds 
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(b) V=24 m-s 
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图 4 不 同 风速 下 路 基 周 围 模 拟 风 沙 流 场 分 布 图 


Fig. 4 Distribution of simulated aeolian sand flow field around the roadbed at different wind speeds 
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于 气流 低速 区 。 而 当 存 在 流动 沙丘 时 ,由 于 流动 沙 
撕 的 阻 潍 耗 能 ,使 得 路 基 坡 面 和 背风 坡 周 围 全 部 处 
于 气流 低速 区 的 范围 ,流动 沙丘 坡 项 处 于 气流 加 速 
区 的 范围 ,此 外 ,路 基 周 边 的 减速 区 范围 明显 增 大 。 
上 述 两 种 情况 下 ,不 同 风 速 的 流 场 分 布 差异 基 
本 相同 , 主要 表现 在 速度 分 区 面积 大 小 的 不 同 ,与 
V=24 m-s }H HE , V=12 ms ! 时 ,其 周边 的 气流 减速 
区 气流 低速 区 面积 减 小 ,气流 加 速 区 面积 增 大 。 
2.2 水 平 风速 分 布 特征 
风沙 流 一 般 贴近 于 地 表 运 动 , 近 地表 风速 对 沙 
粒 运 动 影响 显著 ,研究 近 地 表 沙丘 与 路 基 周 围 气流 
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速度 变化 规律 对 沙 害 的 形成 机 理 具 有 重要 意义 2 。 
选取 不 同 风 速 下 距 地 表 不 同 高 度 (0.3 m、2.3 m、5.3 m, 
6.0m) 处 的 水 平 速度 进行 分 析 ,速度 变化 曲线 如 网 5 
和 图 6 所 示 。 

从 图 5 中 可 以 看 出 , 当 流动 沙丘 不 存在 时 , 距 地 
面 高 度 2.3 mm 处 , 受 路 肩 集 流 效应 的 影响 ,速度 突 
增 ,而 在 路 基 坡 面 ,速度 大 幅 衰 减 ,在 坡 面 中 间 位 
置 ,速度 降 到 最 低 。 对 比 图 5 和 图 6, 当 存在 流动 沙 
ELAN ,沙丘 高 度 以 下 ,气流 速度 在 沙丘 前 方 $~10 m 
范围 内 急剧 减 小 , 遇 到 障碍 物 后 迅速 降 为 0, 远 离 障 
碍 物 影响 后 ,气流 速度 开始 增加 。 在 路 基 背 风 侧 路 
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注 :y 表 示 距 地 表 高 度 。 下 同 。 
图 5 不 同 高 度 下 路 基 周 围 水 平 速度 分 布 图 


Fig. 5 Horizontal velocity distribution around the roadbed at different heights 
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图 6 不 同 高 度 下 流动 沙丘 与 路 基 周 围 水 平 速 度 分 布 


Fig.6 Horizontal velocity distribution of mobile dunes and surrounding roadbed at different heights 
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肩 距 离 地 面 高 度 2.3 m 人 处 ,路 肩 集 流 效 应 的 影响 降 
低 ,气流 增幅 降低 。 沙 丘 高 度 以 上 ,沙丘 顶部 气流 
速度 急剧 增 大 , 受 障碍 物 上 自身 形状 的 影响 ,风速 变 
化 由 “0Q” 型 分 布 变 为 “M” 型 分 布 。 

从 图 7 中 可 以 看 出 ,在 不 同 风速 相同 高 度 处 , 气 
流 的 速度 变化 规律 基本 一 致 。 风 速 分 别 为 12 mes, 
18 ms 24m.s 30m:s 时 ,距离 地 面 $5.3 m 和 6.0m 
处 的 最 大 风速 分 别 为 13.1 m+s',20.2 m-s',27.0 mss", 
34.8 m-s FI 13.3 ms 20.4 ms 27.3 ms 35.2 
mes ,分 别 为 初始 风速 的 1.09、1.12、.1.13、1.16 倍 和 
1.11、1.13、1.14、1.17 倍 , 随 着 风速 增 大 ,沙丘 高 度 以 
上 的 最 大 风速 逐渐 增 大 ,但 幅度 变化 相对 较 小 。 


3 讨论 


当 风 速达 到 临界 起 动 风速 ,在 风力 作用 下 ,地 
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表 沙 粒 开始 起 动 ,脱离 地 表 进 入 气流 中 ,产生 风沙 
运动 ,形成 风沙 流 。 风 沙 流 是 风沙 活动 的 重要 参 
数 ,也 是 造成 风沙 灾害 的 重要 原因 。 随 着 风速 的 增 
加 ,气流 携 沙 量 随 之 增加 ,最 大 输 沙 量 上 升 ,相反 ， 
则 最 大 输 沙 量 下 降 。 图 8 和 图 9 为 不 同 风速 下 流动 
沙丘 和 路 基 周 围 沙 粒 体积 分 数 分 布 图 ,从 图 8 可 以 
看 出 ,流动 沙丘 不 存在 时 ,路 基 迎 风 坡 的 积 沙 量 明显 
大 于 育 风 坡 , 且 随 来 流风 速 的 增加 ,迎风 坡 积 沙 量 逐 
渐 减 少 ,背风 坡 积 沙 量变 化 不 大 ,这 与 石 龙 等 “、 张 
军 平等 ”人 的 人 研究 结 采 相 吻 合 。 从 图 9 可 以 看 出 ， 
存在 流动 沙丘 时 , 积 沙 主要 发 生 在 沙丘 迎风 坡 坡 脚 
处 ,沙丘 背风 坡 及 路 基 周 围 积 沙 较 少 ; 随 着 风速 的 
增加 ,沙丘 迎风 坡 坡 脚 处 积 沙 量 不 断 减少 ,沙丘 背 
风 坡 及 路 基 周 围 积 沙 有 增多 的 趋势 。 对 比 图 8 和 图 
9, 存 在 流动 沙丘 时 , 沙 粒 大 部 分 堆积 在 流动 沙丘 迎 
风 坡 ,路 基 周 围 积 沙 较 少 , 随 着 时 间 的 推移 ,流动 沙 
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图 7 不 同 风速 下 流动 沙丘 与 路 基 周 围 水 平 风 速 分 布 图 


Fig.7 Horizontal wind speed distribution diagram of moving dune and surrounding roadbed under different wind speeds 


(a) V=12 m-s! 


(b) V=24 m-s 


图 8 不 同 风速 下 路 基 沙 粒 体 积分 数 分 布 图 


Fig. 8 Distribution diagram of subgrade sand volume fraction under different wind speeds 


(a) V=12 m-s” 


(b) V=24 m-s 


图 9 不 同 风速 下 流动 沙丘 与 路 基 沙 粒 体 积分 数 分 布 图 


Fig. 9 Distribution of volume fraction of sand grains in mobile dunes and roadbed at different wind speeds 
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丘 附近 的 积 沙 越 来 越 多 ,一 方面 沙丘 本 身 作 为 沙 
源 , 在 沙丘 迎风 坡 会 卷 走 沙 撕 表面 的 部 分 沙 粒 , 并 
沿 坡 面 息 升 , 落 入 沙丘 与 路 基 交 界 处 及 路 面 上 ,为 
一 方面 会 导致 沙丘 慢 慢 沿 主 风 问 移动 ,对 路 基 形 成 
掩埋 ,形成 沙 害 。 

通过 将 数值 模拟 结果 与 图 1 所 示 公 路 路 基 现 场 
只 小 的 对 比 ,模拟 结果 与 现场 实际 积 沙 分 布 基本 一 
致 ,证 明了 数值 模拟 结果 的 准确 性 。 流 动 沙丘 是 路 
基 防 护 过 程 中 不 可 忽视 的 因素 之 一 ,虽然 对 于 穿越 
流动 沙丘 区 的 公路 而 言 ,其 两 侧 的 流动 沙丘 都 进行 
了 表面 沙 障 的 固定 ,但 随 着 时 间 的 推移 , 防 沙 体系 
m TEH, 沙 粒 大 量 在 沙丘 附近 沉积 ,流动 沙丘 会 
慢 慢 向 路 基 靠 近 ,在 路 面 形 成 沙 害 。 因 此 ,对 于 经 
过 流动 沙丘 区 的 公路 , 除 采 取 和 常规 的 草 方 格 BED 
栅栏 和 挡 沙 墙 等 防护 措施 之 外 ,还 应 详细 考察 当地 
的 地 貌 形 态 及 风沙 流 运 动 规律 ,对 沙丘 附近 的 积 沙 
定期 进行 清理 ,及 时 检查 更 新 原 有 的 防护 体系 , 既 
要 阻止 风沙 流 对 沙丘 进行 沙 源 供 给 ,又 要 防止 沙丘 


4 结论 


本 文 利 用 Computational Fluid Dynamics 数值 模 
拟 , 对 不 同 风 速 下 流动 沙丘 区 公路 路 基 的 风沙 流 场 
变化 规律 和 积 沙特 征 进行 了 分 析 , 得 到 如 下 结论 : 

(1) 气流 在 沙丘 和 路 基 附 近 产 生 速 度 分 区 , 当 
存在 流动 沙丘 时 ,路 基 周 边 的 减速 区 范围 明显 增 
大 ,路基 坡 面 和 背风 坡 全 部 处 于 气流 低速 区 的 苑 
围 ,气流 加 速 区 转移 到 流动 沙丘 坡 顶 处 。 

(2) 存在 流动 沙丘 时 ,气流 速度 在 沙丘 顶部 急 
剧 增 大 ,风速 变化 由 "97 型 分 布 变 为 "M 型 分 布 ;不 
同 风 速 下 相同 高 度 处 ,气流 经 过 沙丘 和 路 基 周 转 
时 ,速度 变化 规律 基本 一 致 ,沙丘 高 度 以 上 的 最 大 
风速 逐渐 增 大 ,但 幅度 变化 相对 较 小 。 

(3) 存在 流动 沙丘 时 , 沙 粒 主要 堆积 在 流动 沙 
丘 迎风 坡 ,一 方面 沙丘 本 号 作为 沙 源 ,在 沙丘 迎风 
坡 会 卷 走 沙丘 表面 的 部 分 沙 粒 ,散落 在 路 面 上 , 另 
一 方面 会 导致 沙丘 慢 慢 治 主 风 向 移动 ,对 路 基 形 成 
掩埋 ,形成 沙 害 。 

(4) 位 于 流动 沙丘 区 的 公路 ,应 详细 考察 当地 
的 地 貌 形态 及 风沙 流 运动 规律 ,定期 检查 维护 原 有 


防护 体系 ,及 时 清理 积 沙 , 既 要 阻止 风沙 流 对 沙丘 
进行 沙 源 供给 ,又 要 防止 沙丘 本 号 成 为 沙 源 。 
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Abstract: This study analyzed the wind-sand flow field of a highway roadbed to reveal the formation mechanism 
of pavement sand damage, which provides a theoretical basis for the prevention and control of sand damage on a 
road surface in a shifting dune area. The study area was Highway 303 in Zhongzaohuo, Qinghai Province, China, 
which was selected based on computation fluid dynamic numerical simulation. The numerical simulation results 
showed that when the airflow passes through the shifting dune and the roadbed, the airflow velocity is divided 
into zones. In the presence of a flowing dune, the range of the deceleration zone around the subgrade increased 
significantly, while the subgrade slope and the leeward slope were all in the low speed airflow zone. On the top of 
the sand dunes, the air velocity increased sharply and was affected by the shape of the obstacles, which caused 
the wind speed to change from the “Q” type to “M” type distribution. When the original sand control system on 
the dune surface is not maintained and updated, sand grains accumulate in large quantities near the dune. In this 
manner, the dunes will slowly move in the main wind direction and become the sand source, thus burying the 
roadbed and causing sand damage. The numerical simulation results are consistent with the actual results 
obtained, which indicates the accuracy of the numerical simulation. For the prevention and control of road sand 
damage in the shifting dune area, the local geomorphologic form and the law of wind-sand flow movement 
should be investigated in detail. Checking and maintaining the original protection system regularly, and cleaning 
up the accumulated sand in a timely manner will prevent the wind-sand flow from supplying the sand source, and 
will prevent the sand itself from becoming the sand source. 
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